## ____________________________________________________________________________ ##
##                            PACIFIC HAKE CONTROLS
## ___________________CONTROLS FOR ESTIMATED PARAMETERS________________________ ##
##  Prior descriptions:
##                      -0 uniform (0,0)
##                      -1 normal (p1=mu,p2=sig)
##                      -2 lognormal (p1=log(mu),p2=sig)
##                      -3 beta (p1=alpha,p2=beta)
##                      -4 gamma(p1=alpha,p2=beta)
## ____________________________________________________________________________ ##
7 ## npar

#Original iscam biological oriented priors
## ____________________________________________________________________________ ##

##  ival        lb      ub      phz     prior    p1        p2       parameter name
## ____________________________________________________________________________ ##
     1.9      -1.0     4.0      -1        0     -1.0       4.0      #log_ro  priors - see SS_ro_prior.xls
     0.77      0.2     1.0       1        3      9.766627  2.803034 #steepness  a and b parameters approximate prior from SS - see Betapars.r
    -1.609438 -5.0     0.0      -2        1     -1.609438  0.2      #log.m
     1.9      -5.0    15.0       1        0     -5.0      15.0      #log_avgrec
     1.9      -5.0    15.0       1        0     -5.0      15.0      #log_recinit
     0.15      0.01    0.999     4        3      3.0      12.0      #rho
     1.25      0.01  150.0       3        4      7.49836   5.78354  #varphi (precision)  (RF Change - SJDM had called this kappa)
## ____________________________________________________________________________ ##

#2012 Management oriented priors
##  ival        lb      ub      phz     prior    p1        p2      parameter name
## ____________________________________________________________________________ ##
#   0.2         0.01    3.00      1        2    -1.609438  0.5      #msy -1.609438  0.133939  
#   0.35        0.01    3.00      1        2    -1.049822  0.4      #fmsy
#  -1.481141   -5.0     0.0       2        1    -1.609438  0.1      #log.m   #-1.481141
#   1.163151   -5.0    15.0       1        0    -5.0      15.0      #log_avgrec
#   1.163151   -5.0    15.0       1        0    -5.0      15.0      #log_recinit
#   0.2         0.001   0.999     3        3    12.0      52.8      #rho
#   1.25        0.01   10.0       3        4    39.0625   62.5      #varphi (precision)  

#   0.15        0.01    0.999     4        3     3.0      12.0      #rho
#   1.25        0.01  150.0       3        4     7.49836   5.78354  #varphi (precision)  (RF Change - SJDM had called this kappa)

#   0.223412    0.05    0.9      -1        2    -1.609438  0.1      #m


## _________________________SELECTIVITY PARAMETERS_____________________________ ##
## OPTIONS FOR SELECTIVITY:
##      1) logistic selectivity parameters
##      2) selectivity coefficients
##      3) a constant cubic spline with age-nodes
##      4) a time varying cubic spline with age-nodes
##      5) a time varying bicubic spline with age & year nodes.
##      6) fixed logistic (set isel_type=1, and estimation phase is set to -1 in tpl (ie estimation phase below is ignored))
## Gear 1 fishery:  Gear 2 survey
## isel_type
1 1
## Age at 50% selectivity (logistic) ahat
#4 4.5
#4.82102 4.5
3.5 3.5
## STD at 50% selectivity (logistic) ghat
#1.1 0.5
#1.31762 2.1
0.45 0.45
## No. of age nodes for each gear (0 to ignore).
3 0
## No. of year nodes for each gear (0 to ignore).
5 0
## Estimation phase - any negative number means it is fixed!
1 1
## Penalty weight for 2nd differences w=1/(2*sig^2)
150.0 200.0
## Penalty weight for dome-shaped selectivity 1=1/(2*sig^2)
50.0 200.0
#GAMMA prior for  STD at 50% selectivity (logistic) SAME FOR ALL GEARS - but can switch it on or off
#prior type (4=gamma) par1 par2
#ghat_p1
2.0 4.0
#ghat_pswitch switches the ghat prior on or off - ONE FOR EACH GEAR
0 0
## ____________________________________________________________________________ ##
## ____________________________________________________________________________ ##
##                             Priors for Survey q                              ##
## ____________________________________________________________________________ ##
## nits  #number of surveys
1
## priors 0=uniform density, 1=normal density
1
## prior log(mean)
0.0
## prior sd
0.1
## ____________________________________________________________________________ ##

## _______________________OTHER MISCELLANEOUS CONTROLS cntrl_________________________ ##
0           ## 1  verbose ADMB output (0=off, 1=on)
1           ## 2  recruitment model (1=beverton-holt, 2=ricker)
0.05        ## 3  std in observed catches in first phase.
0.01        ## 4  std in observed catches in last phase.
0           ## 5  Assume unfished in first year (0=FALSE, 1=TRUE)
0.00        ## 6  Minimum proportion to consider in age-proportions for dmvlogistic
0.2         ## 7  Mean fishing mortality for regularizing the estimates of Ft
0.05        ## 8  std in mean fishing mortality in first phase
2.00        ## 9  std in mean fishing mortality in last phase
-1          ## 10 phase for estimating m_deviations (use -1 to turn off mdevs)
0.1         ## 11 std in deviations for natural mortality
12          ## 12 number of estimated nodes for deviations in natural mortality
0.00        ## 13 fraction of total mortality that takes place prior to spawning
1           ## 14 switch for age-composition likelihood (1=dmvlogistic,2=dmultinom)
0           ## 15 1=estimate Management parameters, 0=estimate population parameters
## ____________________________________________________________________________ ##

## RF ADDED NUMBER OF PROJECTION YEARS
##pyrs
3

## RF ADDED harvest control rule switch
##hcr
1 ## 1 = 40-10 Rule ... nothing else implemented yet

##Catch stream from SS (OY)
#SSstream
0.274024
0.282668
0.283845

## eofc
999
